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This talk discusses the power behaved corrections to fragmentation functions for the current 
jet in non-singlet deep inelastic scattering. These corrections are estimated by means of a 
renormalon model using a dispersive approachcl. The assumption of a universal infrared finite 
coupling enables us to provide quantitative estimates. 



1 Introduction 

The experimental study of final state, single hadron, momentum distributions is now being 
carried out at HERA0. Like the energy distribution in e'^e , this quantity naturally depends on 
the parton to hadron fragmentation functions. While it is not possible to calculate fragmentation 
functions within perturbation theory, one can study their evolution with (the logarithmic 
scaling violation) perturbatively and this enables us to extract from the data. However, 
reliance on a purely perturbative approach still causes problems as it does not account for 
higher twist pieces. Though these are suppressed by powers of a large scale, their functional 
dependence on the phase space parameters generally leads to sizable effects, even comparable 
in some cases to the size of the NLO term of the corresponding perturbative estimate. This 
would naturally mean that one must account for these pieces before any comparison with data 
is made. 

The most reliable way in which one can model these effects is the renormalon approach. 
Several papers have been written on these techniques and their applications to structure func- 
tionsi'i fragmentation functions in e'''e~i, and event shape variablesHi. The phenomenological 
success has been encouraging enough to continue applying these techniques. 

For details of the renormalon technique we can refer to the review talks in these proceedings 
0. For our purposes we provide a very brief introduction to the dispersive approach in section 



2. In section 3 we define the observables considered and explain our notation. The calculations 
are described in section 4 and finally quantitative results are presented in the form of a figure. 



2 The Dispersive Approach 

Our estimate of the leading power corrections to the perturbative results are based on the 
approach of reference!. Non-perturbative effects at long distances are assumed to give rise to a 
modification 6aeff{fj?) in the QCD effective coupling at low values of the scale /U^. The effect 
on some observable F is then given by a characteristic function J^{x, e) as follows: 

5Fix, = r ^Saeffif^'mx, e = ^VQ') (1) 
Jo 

where 

j-(x,e) = -e^.F(x,e) (2) 

The characteristic function is obtained by computing the relevant one-loop graphs with a non- 
zero gluon mass fi. 

Arbitrary finite modifications of the effective coupling at low scales would generally introduce 
power corrections of the form l/k'^^ into the ultraviolet behaviour of the running coupling itself. 
Such a modification would destroy the basis of the operator product expansion. This leads to 
the constraint that only terms in the small-e behaviour of the characteristic function that are 
non-analytic at e = will lead to power behaved non-perturbative contributions. In fact we can 
choose to express our results only in terms of without the need to explicitly introduce Oe/j^- 
Then the results can be expressed in terms of ordinary rather than the logarithmic moment 
integrals of the coupling. The relevant moment integral in our case (or more generally for any 
7J2- correction) is 



Cf f°° dfi'^ 
2^ Jo 



M = i^ ^ fj.' Sasifi') . (3) 



3 Fragmentation in DIS 



For the usual reasons (absence of complications due to the remnant jet, similarity with e~^e~) one 
chooses to work only in the current hemisphere of the Breit frame. For more details regarding 
the Breit frame the reader is referred toi and references therein. Since we wish to include only 
particles in the current hemisphere, we define the fragmentation function for a given hadron 
species as a function of the variable z = 2ph.q/Q'^, which measures the fraction of the hadron's 
momentum along the current direction and takes values < z < 1 in the current hemisphere. 
The observable we wish to study is then given by 

where x = Q^/2{P.q) if P is the incoming hadron momentum. The denominator of this expres- 
sion is the fully inclusive deep-inelastic cross section 

5^ = ^{[i + (1-.)Vh.)+2(1-,)F,W} (5) 

where Ft{x) = 2Fi{x) and Fl{x) = F2{x)/x — 2Fi(x) are the usual transverse and longitudinal 
structure functions. The numerator is given by 



dxdQ^dz 



{ [l + (1 - F^{x, z) + 2(1 - y)Fi{x, z)] (6) 



where F^{x,z) = 2F^{x,z) and F^{x,z) = F^{x,z)/x - 2F/* (j;, z) are generalized transverse 
and longitudinal structure functions. The parton model yields F^{x,z) = Fl{x) = and 



F^x, z)=Y: el[q{x)D'^{z) + q{x)D^,{z)] , (7) 
<? 

while 

F!j^{x)=Y.el[q{x)+q{x)]=f{x) (8) 

where and are the quark and antiquark fragmentation functions and q{x),q(x) are the 
corresponding parton distribution functions. The 0{as) result for the numerator in (^) can be 
expressed in terms of 



Ft{x,z) = E.elJ^dU'dC 

x{Ki,ggiC, C)qix/OD^iz/C) + K,,,git Oq{x/OD^Jz/0 



+ Ki,q,{i, C)qix/OD^iz/C) + K,,,giC, Oq{x/i)D^g{z/C) 
+ Ki,gg{C, C)gix/ODhz/0 + Ki^ggiC, C)g{x/i)D'i{z/Q)} 



(9) 



where the F-^ denote the longitudinal or transverse pieces that appear in the numerator on the 
RHS of (4). In the Ki^qq for instance, the i stands for the longitudinal or transverse contribution 
according to what is required and the double suffix denotes the contribution from an incoming as 
well as a fragmenting quark/anti-quark. Similarly the other terms have corresponding suffixes 
that represent incoming gluons (photon-gluon fusion) and fragmenting quarks/anti-quarks or 
incoming quarks/anti-quarks and fragmenting gluons. The variable Q denotes the longitudinal 
momentum of the partons while ^ = with p the incoming parton momentum. In terms of 
the squared matrix elements 

K^,,g{^, C) = l^CpCi^qii, ^ - ^C) (10) 
f\ 

Kr,g,{i, C) = 7rTRCi,g{i, C " ^C) • (H) 

while 

K^,gg{^, C) = ^CpCi^qiC, 1 " C + CO (12) 

where the Ci represent the matrix elements squared! for scattering off incoming quarks or gluons 
as denoted by the suffixes((7 or g). It is emphasised that we do not compute the contribution 
with incoming gluons here. It is expected that for the kinematic range (small x values) at HERA 
the above could be an important effect while perhaps it is not too significant for fixed target 
experiments. 



and 



4 Power corrections 



As mentioned earlier we calculate the power corrections based on the dispersive approach. This 
means computing the first order correction (|9|) with a finite gluon mass e. The coefficient of the 
eloge term then corresponds to the coefficient of the 1/Q^ power behaved term. Computing @ 
with a massive gluon matrix element and phase space and then taking double moments with 
respect to x and z we can express the convolution in that equation as a product involving double 
moments of the K functions and single moments of parton distributions and fragmentation 
functions. Then keeping only the leading non analytic terms in e one finds 



iCT,qg{N,M,e) 
ICL,cc{N,M,e) 

ICL,gg{N,M,e) 

where 



25i (AT) + 2Si(M + 1) - 3 + i + ^ + ^ + ^ 
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-ASliN + 1) - 45i(M + 1) + 6 + 27V + 2M + 2NM 
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(13) 



iV — i -1 

S,{N)=Y.-. (14) 

j=i 3 

and 

£,,,,(iV, Af , e) = /' di [' dC eC^'iC^mi^, (, e) (15) 
Jo Jo 

The /C functions differ from the corresponding K functions we had in (P) only by a factor 
which we choose to include in the definition of our_phenomenological parameter A2. Here we 
have included the virtual contribution calculated inB. The expression given above for the gluon 
fragmentation contribution KLx^qg is valid only for M > 2. There is an infrared divergence 
at M = 0, because we integrate the real gluon contribution over one hemisphere only, which 
does not suffice to cancel the divergent virtual contribution at C = 0. For M = 1 there is a 
contribution of S-^e instead of an e In e term, implying a 1 /Q-correction to this moment of the 
gluon fragmentation function, as is the case in e^e~ annihilation □. For M = 2 the 1 becomes -1 
in the coefficient of elne. All of these changes represent extra contributions at the point C = 0; 
which we can ignore because the fragmentation function at any finite z depends only on the 
behaviour at C > 2; > 0. 

The In e terms generate the logarithmic scaling violations in the structure and fragmentation 
functions while the elne terms give rise to l/Q^ power corrections, as mentioned earlier. We 
also need to consider power corrections in the denominator of and the results for this are 
taken from an earlier publication Bl. 



5 Results and conclusions 



We assume in what follows that the quark fragmentation function is independent of flavour and 
that Dq = Dq which should be reasonable if one confines the discussion to the light quarks and 
if we sum over the fragmentation into all charged particles. One can combine the leading order 
results ^ with the power corrections and express the results as 

F\z-x,Q^) = Dq{z) + ^j^SldUldU{x/0{[HT,,,{i.C) , . 

where f{x) is the charge- weighted parton distribution and HT^q{C) is the higher- twist coef- 
ficient function for the transverse structure function 

Ht,,{0 = ^^^^ - 2 - 4e + 45(1 - + - ■ (17) 

"Note that the definition here difTers from that in Refs.E by a factor of — 1/^. 




The negative contribution from the structure function comes about because of our normalization 
to the structure functions in equation Also the longitudinal contribution is identical between 
the numerator and denominator of equation which means it does not appear in the final result. 
We represent our final result as 

F\z;x,Q^) = Dgiz) + 

We plot the function H{z, x) as a function of z for different values of x. 

For the plot we use the ALEPH0 parametrizations of the light quark and gluon fragmen- 
tation functions for charged hadrons at Q = 22 GeV, and the corresponding MRST (central 
gluon) EJ parton distributions. Thus the predictions are at = 484 GeV^, but H{z; x) depends 
only weakly (logarithmically) on Q^, and in any case our method is not reliable at the level of 
logarithmic variations. Results become insensitive to x below the values shown in Fig.l. Recall, 
however, that we have not computed the singlet contribution, which may well be important at 
low X because of the increase in the gluon distribution there. 

The predicted power corrections are qualitatively similar to those for fragmentation functions 
in e~^e~ annihilation Q, though somewhat larger in magnitude. Part of the increase comes from 
the negative higher-twist correction to the transverse structure function in the denominatoii of 
Eq.(^). The contribution from gluon fragmentation, although subject to further correctionsQ is 
estimated to be relatively small for z > 0.2. 

The last point that needs to be made is that inclusion of electroweak effects does not change 



(18) 



the qualitative situation and in fact just means redefining the function /(x) to include elec- 
troweak couplings rather than just the quark charges. 
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